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ABSTRACT 

It is argued that renormalisation group flow can be interpreted as being a Hamiltonian 
vector flow on a phase space which consists of the couplings of the theory and their conju- 
gate "momenta" , which are the vacuum expectation values of the corresponding composite 
operators. The Hamiltonian is linear in the conjugate variables and can be identifled with 
the vacuum expectation value of the trace of the energy-momentum operator. For theories 
with massive couplings the identity operator plays a central role and its associated cou- 
pling gives rise to a potential in the flow equations. The evolution of any quantity , such 
as A^-point Green functions, under renormalisation group flow can be obtained from its 
Poisson bracket with the Hamiltonian. Ward identities can be represented as constants of 
the motion which act as symmetry generators on the phase space via the Poisson bracket 
structure. 
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§1 Introduction 

The history of relativistic quantum field theory is plagued with divergences. The 
canonical cure for this is to replace the (infinite) bare co-ordinates Qq by (finite) renor- 
malised co-ordintates g%{n) at some renormalisation point k. This requires regulating the 
divergences and choosing a subtraction prescription by introducing counter terms. Even for 
theories with no divergences it is sometimes useful to introduce counter terms and define 
renormalised couplings - an example of this is the e-expansion around the Wilson-Fisher 
non-trivial fixed point in 3-dimensional massive theory [1] 

However knowledge of the renormalised couplings alone is not sufficient, one must also 
know what the counter terms are. This is (almost) equivalent to knowing the /3-functions 
of the theory, since the counter terms are specified by the difference ^g"'{K) = Qq — fi'^(«)- 
Demanding that the bare couplings are independent of the renormalisation point gives 

oa_d9'h_ d{Ag-) 
^ ~ dt ~ dt ' 

where t = lnK. 

The traditional approach has been first to choose the counter terms and then calculate 
the /9-functions. However the counter terms contain a certain ambiguity, they can be 
modified by adding a finite function of the couplings to Ag"'. Thus, as is well known, the 
/3-functions are not unique, they depend on the subtraction procedure. In practice though 
they are constrained by the requirement that a renormalisation prescription be chosen 
which leads to a perturbation theory that makes sense i.e. converges reasonably quickly, 
at least asymptotically. Of course, the /3-functions are not completely arbitrary - they 
have zeros which cannot be removed by changing prescription. Viewing them as vector 
fields ( "velocities" ) on the space of couplings, a zero of the vector field in one prescription 
remains a zero in any other prescription (of course individual components of the vector (3 
may vanish in one prescription and not in another, but a zero of the vector field requires 
the entire vector to vanish). Thus, in one sense at least, a change in renormalisation 
prescription can be thought of as a co-ordinate transformation (diffeomorphism) on the 
space of couplings since this also leaves the zeros of a vector field unchanged (although it 
may, and in general would, change the numerical values of the co-ordinates of the point at 
which the zero occurs). 

One can imagine turning the logic round and first choosing the /^-functions and then 
using equation (1) to determine the counter terms (up to an arbitrary constant i.e. a 
renormalisation group invariant). Of course the choice must be judicious - it would be 
crazy to choose a positive /9-function for massless QCD, the resulting theory would be 
completely unstable and probably would not even exist. In theories with more than one 
coupling, however, there is more freedom. If one were unlucky enough to choose a j3- 
function which pointed away from a nearby attractive fixed point the RG evolution would 
presumably force it to turn round and point in towards the fixed point - behaviour which 
cannot happen if there is only one coupling. Choosing a /3-function with the "wrong sign" 
in a theory with more than one coupling is not pathological, it is merely an indication 
of a "bad" choice of co-ordinates. Transforming to "good" co-ordinates (i.e. co-ordinates 
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which have a sensible physical interpretation and are not just abstract parameters) would 
require performing a diffeomorphism which "unwinds" the RG trajectory so that it flows 
into the flxed point in a straight line. In actual fact, from a topological point of view, the 
notion of "pointing away" does not make sense unless a metric is defined on the space of 
couplings. To have a concept of "pointing away" requires having a definition of angles, 
and this needs a metric. Without a metric any direction is much the same as any other, it 
is difficult to navigate in a space without a metric! 

Another quantity which should not change under a co-ordinate transformation is the 
signature of the matrix of anomalous dimensions, da(3^, at a fixed point. The number of 
positive and negative eigenvalues of this matrix determines in how many directions the RG 
flow is attracted to the flxed point and in how many it is repelled from it. For massless 
QCD this matrix is one dimensional and this is another way to see that the sign of the 
/3-function cannot be changed by a co-ordinate transformation in a theory with only one 
coupling. 

There is an analogy here with classical mechanics. Consider a theory with n — 1 
couplings, g"', a = 1, . . . ,n—l (from now on the subscript R on couplings will be dropped in 
this section - all couplings are renormalised unless otherewise indicated) . Denote the space 
of couplings by Ai. The /3-functions P°'{g) constitute a vector field on the n — 1-dimensional 
differentiable manifold Ai (for the moment the topology and global properties of Ai will not 
be relevant and it may as well be taken to be R*^"-*^ - consideration of the global topology 
will be restricted to a few comments in the final section). The 2n— 2 dimensional space with 
co-ordinates {g°',P°') (the tangent bundle T{M.)) is thus analogous to the configuration 
space of co-ordinates and velocities in classical mechanics. Choose a point in T{Ai) and 
the evolution of the system is determined by the dynamics. All the necessary information 
for computing the RG evolution is contained in the generating functional (or free energy 
in statistical mechanics) W{g,t) = J w{g,t)d^x = — \nZ where w{g,t) is the free energy 
density. * Just as in classical mechanics one can imagine transforming to a phase space, 
with co-ordinates ((7", (f)a) (the cotangent bundle T*{M.)) where are "momenta" dual to 
the "velocities" (3°'. A natural choice for the (f)a is the vacuum expectation value of the (in 
general composite) operator associated with the coupling g"-, (pa — ^^ga'*^ [2] . In order to 
streamline some of the formula it will be convenient to re-scale all the couplings by their 
canonical dimensions so that they are dimensionless. When this is done the variables 0a 
are densities with canonical mass dimension D. 

Following O'Connor and Stephens [2] one can ask if there might be some notion of a 
Hamiltonian function on phase space H^g^cp^t) which could govern the RG evolution of 
the couplings and the expectation values together so that the RG flow can be regarded 
as a Hamiltonian vector flow on phase space. The answer is yes and the construction 
turns out to be remarkably simple. A Hamiltonian will be presented in section three which 
is linear in the momenta, rather than quadratic as in non-relativist ic particle mechanics. 
It is in fact minus the expectation value of the trace of the energy-momentum operator. 



* For simplicity we shall work on flat D-dimensional Euclidean space, so that transla- 
tional invariance ensures that w{g,t) is independent of position, but the concepts presented 
here are more general than this. 
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H = — < T >. Despite the linearity of the Hamiltonian in the momenta, the dynamical 
evolution is not trivial and some examples will be examined to show this in detail. 

The construction requires a careful consideration of the role of the identity operator 
in theories which involve massive couplings. This can be done by introducing a coupling 
associated with the identity operator, A (a cosmological constant). Its conjugate momen- 
tum, (/>A, is just the expectation value of the identity operator. The /3-function associated 
with A turns out to be linear in A [3] , /3^((7, A) = —DA + U{g) where U{g) is an analytic 
function of the other couplings and is independent of A. The conventions adopted in this 
paper are such that the cosmological constant is scaled by its canonical dimensions and 
A is dimensionless. This explains the term —DA here which not otherwise be present. A 
consequence of this is that the 0a = is also a density. For massless theories, U (g) — 0. 
For the purposes of the introductory discussion presented here it will be assumed that the 
/3-functions have no explicit k dependence and only depend on k implicitly through 
A full generalisation to the situation where a subtraction procedure is chosen in which the 
/^-functions have explicit k dependence is given later. 

It will be shown in section three that the Hamiltonian defined by 

H{g,<j>)^P'^{g)<j>a + P\g,A)<l>A, (2) 



governs the renormalisation group flow of the couplings g°' and the expectation values (j)a- 
Of course, once the theory has been solved, all of the vacuum expectation values (VEV's) 
can be expressed as functions of g"' and t, (pa = 4>a{9-,t) and then the Hamiltonian can 
indeed then be written as a function of 5^", A and t alone (for example 0a = n^). However 
the philosophy here is that at the outset g"" and 4>a are to be considered as independent 
variables and H{g,(j)) depends on each separately. 

The main result of this paper is that the RG evolution of g"' and 0^ is given by 
"Hamilton's equations" , 

dg'' dH 



dt dcpa _ 
rf0a dH 



dt dg"" 



(3) 







The first equation follows simply from the definition of H{g, 0) in (2) while the second 
contains non-trivial dynamics, despite the simple form of H. Indeed one can interpret 
U (g) as a potential and re-write the second of equations (3) as 

^ = -K^dUig), (4) 

where = (padg"' is a one-form and dU = ^^dg°' is the exterior derivative of the poten- 
tial. The analogy between (4) and Newton's second law, ^ = —dU for a particle with 
momentum p moving in a potential U, is obvious. 

The idea that RG flow might be related to a potential was flrst suggested, to the au- 
thor's knowledge, by Wallace and Zia [4] but so far the investigations in this direction seem 
to have been attempts to flnd a potential for the /^-functions which requires introducing 
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a metric on the space of couplings, rather than for the VEV's as suggested here. In the 
construction presented here a metric on T{Ai) is not necessary. 

The reformulation of the RG evolution in terms of Hamiltonian flow allows the intro- 
duction of Poisson brackets and their associated symplectic structure. (That a symplectic 
structure should be relevant to RG flow was first suggested, to the authors knowledge, in 
[2].) Extending the set {g"'} to include the cosmological constant, {(7", A} = {g°^} where 
a = 1, ... ,71, the Poisson bracket of any two quantities A and B is, of course, given by 
expressing them as functions of g"' and ^a, possibly also with an explicit t dependence, 
and taking the combination 

dg^d4>a d4>adg^- 



Obviously 



{/, = Sa'. (6) 



The RG evolution for any function on phase space is then given by 

dA_dA 
'dt ~ ~dt 



+ {AH}. (7) 



Note in particular that, when there is no explicit k dependence in the /^-functions, 
the Hamiltonian (2) is a constant of the motion ^ = 0. Once more we stress that, even 
though it is trivial that (f)A = it would be wrong to include this k dependence in if ((7, (j)) 
explicitly - (p\ is to be considered as an independent variable in this formalism and H has 
no explicit k dependence. One could, however, omit A as an independent variable and 
define a t-dependent Hamiltonian h{g,(l),t) = f3"'{g)(j)a + K^U{g) on T*{M). This leads 
to the same equations of motion (3) but h{g, (j), t) is not a constant if there are massive 
couplings in the theory, instead ^ = ^ = Dn^U{g). 

The classical analogy can be taken even further. It is also argued in section three that 
the RG equation for the generating function w{g^ t) can be written as 



dw 
'dt 



^ + h(..^)=0, (8) 



which is clearly a field theoretic version of the Hamilton-Jacobi equation. 

The layout of the paper will be as follows. In §2 the RG equation for w{g, t) is 
derived, taking particular care over the role of the identity operator and the way that 
the cosmological constant is related to masses. In §3 a Hamiltonian formalism of RG 
evolution is developed. A symplectic structure on the phase space {g"',(j)a) is introduced 
and Hamilton's equations (3) are derived together with the Hamilton-Jacobi equation (8). 
The renormalisation group equation for A/"-point Green functions is presented as a special 
case of equation (7). §4 (which is the only part of this paper which uses perturbation 
theory) exemplifies the ideas with the use of massive A(^^ in four dimensions as a model. 
The implementation of symmetries is discussed in §5, where some examples are used to 
argue that Ward identities can be represented by symmetries on phase space and can 
be used to construct RG invariants ("constants of motion") which commute with the 
Hamiltonian and generate the symmetry through the Poisson bracket structure. Finally 
the results are summarised in §6. 
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§2 The Renormalisation Group Equation For The Partition Function 



Consider a renormalisable field theory, parameterised by a set of renormalised cou- 
plings Qj^, a = l,...,n — 1 where n is finite (the superscript n will be reserved for a 
cosmological constant). In this section all renormalised quantities will be denoted by the 
letter R so as to try to keep the argument as clear as possible. The bare couplings Qq can 
be thought of as functions ggigR, e) of the renormalised couplings plus a regular isat ion pa- 
rameter, e (strictly speaking these functions should also contain an explicit k dependence 
since their total derivative with respect to k must vanish, but this is not shown here). The 
bare couplings are analytic functions of the (7^ provided e ^ 0, but are singular in the limit 
e — > 0. The transformation gQ — > gf^ can be viewed as a co-ordinate transformation on the 
n — 1-dimensional space of theories, A4. This co-ordinate transformation is singular in the 
limit e but, as long as the theory is renormalisable, this is not a disaster and can be 
treated consistently. For the moment the method of regularisation is left open, one could 
for example use a cut-off Ac and set e = k/Ac or dimensional continuation with e = D — 4. 

The action can be written as a linear combination of "basic" operators ^»°, which will 
include composite operators, 

So{go,^^e) = J Co{go,^\x))d''x where Co = g^^l{x). (9) 

Strictly speaking, since space is taken to be Euclidean, this is the energy rather than the 
action - but a Hamiltonian will appear in the next section in a totally different context. 
By abuse of language therefore 5*0 will be called the action. 
In massive Xip'^ theory, for example, one would have 

A = kod^ipod^ipo + joipo + ^ml(fl + ^(fg (10) 

with $° = d^ipod^ipo, = (fio, $^2 = and $^ — ^(/Jq and four independent couplings 
koijoiiT^o ^-nd Aq. Without loss of generality the fields can be re-scaled to set A;o = 1- If 
necessary the couplings can be made functions of position. Thus the notation can be 
extended to allow source terms, jQ{x)(fQ{x), with jo{x) a function of position which may 
be set to zero after all (functional) differentiations have been carried out. More generally 
all of the couplings can be made to depend on position so as to introduce sources for the 
composite operators as well [5] . After all differentiations have been carried out these 
sources can be set to constant values if so desired. 

Returning to the general case, consider a renormalisable theory written in terms of the 
bare couplings. Representing all of the bare fields generically by (po, the partition function 
(generating functional) is, 

Zo{go) = J V<poe-^°^^°'^°\ (11) 

The generating functional for connected Green functions (the free energy) is Wo{go) = 
— InZo(^o)- As usual the free energy density is defined via Wo(^o) = / wo{go,x)d^x. The 
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introduction of wq avoids trivial volume divergences in the case of infinite space. It has 
canonical mass dimension D. 

For future convenience a coupling for the identity operator will be included, 

Z{go,Ao) = J P^oe-^°^^°'*°)-/''""^° = e" / '^°^^°Zo(^o) (12) 

(there is no bare subscript on Z{gQ, Aq) because, as will be explained later, it is finite). In 
analogy with general relativity Aq might be called a cosmological constant, but in different 
physical situations it would have different physical interpretations. 

It is important to realise that Aq is independent of Qq and plays the role of a new 
coupling for the identity operator. Thus the set {qq} can be extended to {Qq} = {qqjAq} 
where a = 1, .... n are co-ordinates on a n-dimensional manifold Ai. Aq can be included 
in the bare Lagrangian as Cq — Qq^^ + ^ol- In terms of densities 

W{go,Ao) = -lnZ{go,Ao) = J w{go, Ao)d'^ x, (13) 

with w{go,Ao) = wo{gQ) + Aq linear in Aq. 

When all couplings are independent of position and the theory is translationally in- 
variant w{go, Aq) is independent of x and this will be assumed from now on. In situations 
in which translational invariance is not a symmetry of D-dimensional space, one will need 
to introduce extra terms involving the Riemann tensor into the action [3] . 

Expectation values of bare quantities can be obtained by differentiating w{go,AQ) with 
respect to the couplings, 

< >= ( — — ) = - — and 1 =< 1 >= — — . (14) 
\dgS/ dg^ dAo ^ ' 

The bare operators = |^ are co-vectors on the space of couplings, i.e. ^» := ^^^dg^ 
is an operator valued one-form. This notion of the basic operators being co-vectors on 
the space of couplings is implicit in the work of Zamolodchikov [6] . For a conformal field 
theory in two dimensions would be the primary fields of the theory. Expectation values, 
dw —< $9 > dgQ, are (exact) real valued one-forms on M.. 

The bare operators can be written as linear combinations of renormalised operators 

= Za'^f, (15) 

where Zq,^ is a matrix of renormalisation constants (see for example [7] ). In general there 

will be operator mixing and Za^ will not be diagonal. If the bare action has no massive 
couplings it is not necessary to include the identity operator in the list and a cosmological 
constant can be omitted, since Za^ — {a < n — 1) when all couplings are massless. But if 
there are massive couplings in the bare action then some or all of the Za"' will be non-zero 
and Aq plays a crucial role. 
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Following up the idea that the transition from bare to renormalised couplings can be 
implemented as a co-ordinate transformation on Ai, we consider gQ^gn^e) to be analytic 
functions of (keeping e for the moment). Consider 

dg%dgl dg% ^" ^ ^ 

Clearly these should be related to the renormalised operators ^f. In fact 

is the inverse of the operator mixing matrix Za = g|f [8] . Thus 

g = (Z-')/4.? = *?. (18) 

The operator valued one-form $ can now be expressed in either co-ordinate system, 

$ = ^Idgt = ^&dg% (19) 
and similarly the real valued one-form dw is 

dw =< $ >=< $5 > dg^ =< $f > (20) 

Thus 

g = ^ < *S >= < >=< *f > ■ (21) 

Partial derivatives here mean, of course, that all quantities are considered to be functions 
of the renormalised couplings, and each is varied independently of the others. When 
massive couplings are present Aq must be considered to be a function of the renormalised 
couplings. This is necessary because the $o then mix with the identity operator under 
renormalisation. 

Note that, since both g^^ and < > in equation (21) are finite as the regularisation 
parameter e — >^ 0, w{gR, t) considered as a function of the renormalised couplings and the 
renormalisation point k must also be finite as e — 0. This is why there is no bare subscript 
on the definition of Z{gQ, Aq) in equation (12), [3], and is one of the reasons for introducing 
a coupling for the identity operator - both Zo(fi'o) -^o seperately diverge as e — > 

but the combination Z{g%) = Z{g^) = e~ I '^'^ ""^'^ Zo{gS) is finite. A finite, renormalised 
generating function, Wfl(5r^,t) = J WR{gj^,t)d^x = — In Zii{gj^), can now be defined by 

Z{9l) = e-/'^"^^°Zo(^o") = Z{g%,t) = e'^" I ZM,t) (22) 
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(the factor in the exponential is in accord with the convention that is dimensionless) . 
Thus w = 'u;o(fi'o) + = wjiigf^, t) + Aj^k^ can now be interpreted as a finite function of 
the renormahsed couphngs with the crucial property that 

=< > and ^ =< $^ >=< >= Av^. (23) 



Note that Aq is linear in Aji, a fact which follows from the observation that the identity 
operator 1 = does not get renormahsed. Thus 

Indeed even the /^-function for Ar only depends on Ar linearly through canonical dimen- 
sions, since writing 

Ao = Av^(A« + F(^^,e)), (25) 

we have 

= ^ = -DAr + U{g-R) (26) 

where U{g'j^ ■~ depends only on the gf^ for a < n — 1, not on Ar, and is finite as 
e — > 0. The fact that U{gR) is independent of A can be seen from the following argument. 
U{gR) is a quantum correction to the canonical dimensions of the cosmological constant 
and as such can be determined (in principle) using perturbation theory and Feynman 
diagrams, but A cancels out of all Feynamn diagrams due to the normalisation factor ^. 
Hence U (gR) is independent of A. 

To make contact with the familiar notions of a perturbative analysis one can write 
the bare couplings in terms of the renormahsed couplings as 

gl{gR,e)^gi + Ag%gR,e) (27) 

where A^r'* is a correction (which diverges as e — > 0), so that 



Co = g^^l = ig% + Ag^)Za'^^ 

a 



(28) 

)f + counter terms (C.T.'s). 
This gives 

= Z(^o, Ao) = e-/'^"^^° I P<^oe-/'"^^«*"(^)+^■^•'^ (29) 



—w 

e 



where the identity operator has been included among the counter terms. Absorbing a 
further term proportional to the identity into the counter terms we have 



—w 

e 



If desired, one can perform the functional integral over renormalised fields rather than 



bare fields by setting VipQ — ^ T^^r^ where lpq 



.1/2 



ipji with z being the wave function 



renormalisation factor, and then adding a further term ^ In^ J d^xl to the counter terms. 

Leaving perturbation theory behind and returning to the general analysis, we are now 
in a position to write the renormalisation group equation for the free energy, w{g'^^t). 
Since w = wo{gQ) + Aq and all bare couplings are independent of the renormalisation point 
t = \nK, we have 



dw{g%,AR,t) 
dt 



dw 

d?R 



R 



dt 



9r,-^r 



If we now denote the VEV's by 0^ :=< >, this reads (since 



dw 



(31) 



K,^ = (j)^ and 



dw 



) 



dw{g%,t) 



dt 



, +/3"(^7S)0f + /3^(^7|)0A = O. (32) 
9r,Ar 

This is the equation that will be used in the next section to argue for Hamiltonian flow on 

(9%, <P&) space. 

Alternatively, since A^^ only ever appears in this equation linearly, it can be eliminated 
using equation (26) and w = WR{gfi) + Ark^ , with wr independent of Ar, to give 



dwR{g%,t) 



dt 



9r 



(33) 



Note that it is not true, in general, that ^^'^^^R'^^ — o, since the presence of massive 
couplings necessitates the introduction the function U{g) in a general renormalisation 
prescription. 

This analysis has been a somewhat lengthy treatment of concepts that are not new, 
but it has been included in order to expose clearly the role of the cosmological constant in 
theories with massive couplings as well as to set up the notation. 
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§3 Symplectic Structure And Hamiltonian Flow Of The RG Equation 



In this section it will be shown that the RG equation derived for the generating func- 
tion in the previous section naturally admits a symplectic structure with its concomitant 
Poisson brackets, and the renormalisation group flow can be obtained from a Hamilto- 
nian function on phase space in a manner analogous to dynamical evolution in classical 
mechanics (but with important differences). 

The starting point is equation (32) 



^+U'=''- '''' 



From now on the qualifier R on g'^ and (j)?^ will be omitted as all quantities will be renor- 
malised, unless otherwise indicated. The generating function w appearing in this section 
is always w = wji + k^Aji. It is stressed that w is linear in A. 

To highlight the analogy with classical mechanics, we shall define a function H, 

so that (34) can be written 

H{g,4>)=(3\g)<t>a + (^''<t>K. (36) 

(For the moment it will be assumed that the /^-functions have no explicit k dependence so 
that H{g,(f)) has no explicit t dependence - a generalisation including explicitly k dependent 
/3-functions is given towards the end of this section.) The philosophy now is to forget where 
the (pa came from and treat them as independent variables. It is only after the theory has 
been solved that we can use (pa — daW. 

Consider the left hand side of equation (34) as a differential 

G = (f)i,dg^ - Hdt, (37) 

where Q{g, </>, t) is a one-form on the 2n + 1 dimensional space parameterised by g"', (pa and 
t. When the theory is solved, and (paig^t) is written as an explicit function of g°' and t, 
= dw is exact, but when the (pa are treated as independent variables is not exact. The 
discussion now parallel's the treatment of classical mechanics in [9] . Just as the couplings 
g"'{t) evolve along the RG trajectories, so do the expectation values (pa{t)- Thus the RG 
trajectories can be pictured as flow lines in (gf", (pa, t) space. Treating (pa as independent 
variables construct the two-form 



9. = -nijdx' A dx^ = dQ = d(pa A dg^ -dH Adt 



1 

_( 

dH 

d(pa A dg^ - 



dg^ 



dg" Adt- 



(p d(pa 



d(pa A dt, 

9 



(38) 
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where {x^} = {g"", (f)a, t} i = 1, . . . , 2n + 1 are co-ordinates on the 2n+l dimensional space. 
Now Q can be written as an anti-symmetric matrix, 




I 



dH_\ 
d(f>a\ 



dH \ 
dH I 
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\ 



I 



(39) 



where I is the n x n identity matrix. Since is an odd dimensional anti-symmetric matrix 
it must have at least one zero eigenvalue (it will be assumed that it has only one, otherwise 
the restriction of O to surfaces of constant t would result in a degenerate symplectic form 
on phase space). The corresponding eigenvector, ^, is easily seen to be 



dH 



99" 



1 



(40) 



It seems natural to identify the flow lines of the vector field ^ with renormalisation group 
trajectories, since © is exact when the theory is solved and (j)a{g,t) are substituted into 
equation (37). This requires 



dg'^ _ dH 

dt d(f)a 



and 



OH 



9 



dt 



(41) 



Obviously ^ = /3" by definition, so the first equation is certainly consistent. For d 



n. 



the second equation reduces to the identity 



Dk 



D 



— . The interpretation of the other 
equation hinges on the crucial observation that the function U{g) — P^ + DA is an analytic 
function of the g"^ (independent of A). Thus, from equation (36) 



d(f>a 
dt 



d 



dg 



-a{-''U{g))- 



'dg^ 



or 



dt 



(42) 



(43) 



This is the renormalisation group equation for the RG evolution of the vacuum expectation 
values of the basic operators of the theory. The parallel with Newton's second law is 
obvious. The matrix of anomalous dimensions dafi'^ appears as a pseudo-force (Coriolis 
force) and the function U {g) is a potential. For massless theories U vanishes, thus massless 
theories are analogous to free particle motion. Just as in classical mechanics the Coriolis 
force can be eliminated when the motion of the basis vectors is included in the equation. 
Consider therefore the one-form = (j)adg°'. One has 



dt 



, a 1 didg"") 
dtj''^^^^ 

+ hTr:;]dg'', 



dt 



dg'^ 



(44) 
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where we have used ^^^-^ = d{^) = dhl3°'dg^. Thus equation (43) can be written in 
co-ordinate free notation as, 

-K^'dUig), (4) 



dt 



where dU = ^^dg"^ is the exterior derivative. Of course, once the theory is solved, the (pa 
can be expressed as explicit functions of g"' and t, (f)a{g,t), so that 



d0 
'dt 



and equation (43) then becomes 

d<P. 



dt 



+ 



p'^dbcl^a + {da(3')(l)b] dg\ 



(45) 



dt 



9 



(46) 



which is a version of the RG equation for the VEV's, including the anomalous dimensions 
and the inhomogeneous term —K^dall which arises due to masses. 

Yet another way of expressing this is to observe that the left hand side of (46) involves 
the Lie derivative, C0(j), of the one- form (padg"' with respect to the vector field (3 = (3^-^ 
[10] . Since is exact we have 



C^^ = d{i0cf>) = —^^(3'ig)M9,t)}dg-, 



(47) 



where i^cj) denotes the contraction of /5 with the one- form 0, = /?"0a- Thus another 
way of writing (43) is 

= -d{K^U{g)+p\g)M9)}, (48) 

where again d = dg"--^. 

The analogy with classical mechanics can be taken further still. The definition of H 
in equation (35), 

^(^,0) + ^ = O, (49) 

can be expressed as a partial differential equation in the n + 1 variables (^r", t). Since, when 
the theory is solved, 4>a = we have 



dw 
'dt 



9 



( dw \ 



(8) 



which is clearly an analogue of the Hamilton- Jacobi equation. (That the RG equation 
ought to be expressable as a Hamilton- Jacobi type equation was first suggested to the 
author by Denjoe O'Connor and Chris Stephens [11] .) Thus the generating functional in 
quantum field theory (or free energy density in statistical mechanics) is playing the role of 
the action in classical mechanics (Hamilton's principal function). 
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This structure suggests a reformulation of the renormahsation group. Instead of ex- 
pressing the RG running in terms of /3-functions and couphngs, which can be thought of 
as co-ordinates on configuration space (the tangent bundle T{Ai)) it may be useful to use 
instead phase space variables (the co-tangent bundle T*{/A)). Any quantity, A, should 
then be considered to be a function of the 2n co-ordinates ((7", (pa) and possibly also the 
renormahsation point t = \nK. The RG evolution of A{g, </>, t) is then given by, 

dA _ dA dH 
Itt ~ WWa 



dA dH dA 



9,4> 



(50) 



where {A, H} is the usual Poisson bracket with {gr", 0^} = 5"^. Since there is no explicit k, 
dependence in the Hamiltonian (36) H{g, cj)) is a RG invariant (a constant of the motion) 



dH ^ 



(61) 



but this is only true when there is no explicit k dependence in the /^-functions. 

In particular the RG evolution of iV-point Green functions is of key importance in any 
theory. These can be viewed as rank N tensors on the space of couplings, 

Ci^l.aj^ {Xi,..., Xn) =< ^ai (Xi) • • • ^ajv (xn) >, (52) 

where ^ai{xi) = (x^) — (pai has zero vacuum expectation value (and is independent of 
A). The RG equation for ci^l.a^ {xi, . . . ,xn) is obtained by the process described above, 

d 



dt 



GiZ-aN {xi,...,Xn) 



d_ 
dt 



(P.t 



d 



d/3' 

d?J d(Pb 



g, t d(ph 



9,t 



G 



{N) 



{xi,...,xn)- 

(53) 

If the co-vector basis dg"' is also included, so as to write the tensor in co-ordinate free 
notation G^^'' — Gi^^.ajv(^i; • • • ; XN)dg°'^ ■ ■ ■ dg"'^ , one arrives at the equation 



rdG^\ 

V dt J 



- ^G 



d 



+ J2(9a.f3')G 

ai---ai-i 



N 



dgb 



cP,t 



(N) 



ba. 



+i---«jv 



i=l 



d 



dcPt 



G 



(N) 



9,t 



0,1 •■■aiv 



d 



dcPi 



9it 



(54) 

The RG equation for A/"-point Green functions in the Hamiltonian formalism is finally 
obtained by observing that G^^^ can equally well be written in bare co-ordinates and so 

is independent of t. The left hand side of (54) therefore vanishes, = 0, and the RG 

equation is 



^Gif^la^ {g, 4>, t) + (3%GZla^ {g, 0, t) + $^(a„,/3'')Gif.?.,^_^,„^^^...,^ {g, 0, t) 

i=l 

= {MdhP"") + i^^'dbU) Q^Gi'ila^ {g, <P, t). 



N 



(55) 
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Alternatively, using equation (46) and allowing for the g°' and t dependence of 0a after the 
theory has been solved, this equation can be re-expressed as 

Were it not for the last term on the left hand side of this equation, it would just be the 
definition of the Lie derivative of G^^^ with respect to the vector field (3 - the last term 
is a correction to this interpretation. That the RG equation for A'"-point amplitudes could 
be written as a Lie derivative was observed in [10], and corrections to this interpretation 
were investigated in [12] . 

Returning now to the general formalism let us consider more general canonical trans- 
formations. A renormalisation group transformation is like time evolution in classical 
mechanics and as such preserves the symplectic structure on the 2n-dimensional space 

^ = ^lt=const = d(l)a A dg". (57) 

Of course, there will in general be other important canonical transformations which do 
not necessarily correspond to RG transformations. In particular, the transformation from 
bare to renormalised couplings preserves cD, 

UJ = d(|)!,^dg^ = d4>l^dgl (58) 

but in a general renormalisation prescription this would not correspond to a RG transfor- 
mation (except perhaps if the theory is regularised by using a cut-off and BPHZ renor- 
malisation, in which case the bare couplings do have the interpretation of just being the 
renormalised couplings at some very high energy, Ac). Since the bare couplings are in- 
dependent of the renormalisation point, = = 0, this transformation is analogous 
to the canonical transformation in classical mechanics which takes one from time depen- 
dent phase space variables {q°' (t) , Pa{t)) to the initial point (q'Q,7>°). Referring back to the 
Hamilton-Jacobi equation (8) wc see that the generating function for this canonical trans- 
formation, in the familiar classical mechanical sense, is nothing other than the generating 
function(al) of the quantum field theory, w. 

The analogy with the Hamilton-Jacobi equation of classical mechanics can be further 
highlighted by explicitly indicating that w{g, t) depends on the subtraction procedure and 
writing it as w{g{t), go, t). This emphasises its dependence on the counterterms Ag"- = g^ — 
g°'{t), and does not affect the argument that ^ = 0. In BPHZ renormalisation gg really 
can be thought as lying on on the RG trajectory and the analogy between w{g(t), go,t) 
and Hamilton's principal function, the action S{q, qo, t) along a classical trajectory is even 
stronger. 

A crucial difference between the phase space approach to the RG presented here and 
classical mechanics lies in the Legendre transform, 

iy(^,0)-/3>a = O, (59) 
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with = ■§^\g, which vanishes and in particular is not invertible. However, as the 
examples of the next two sections show, the flow is still far from trivial! 

The Legendre transform presented here is also much simpler in form than quantum 
field theory Legendre transform introduced by Jona-Lasinio [13] .* The latter is imple- 
mented at the level of the generating functional w itself, rather than on its derivative 
and this leads to the effective action, which most certainly does not vanish. 

Note that the Hamiltonian (36), 

= + /?^0A, (60) 

actually has a simple physical interpretation. The right hand side of this equation is just 
the negative of the usual definition of the vacuum expectation value of the trace of the 
energy-momentum tensor of the theory, H = — < T >. It should not come as a surprise 
that < T >= ^1 g since varying t with the couplings fixed is completely equivalent to a 
conformal rescaling of the metric. The derivative ^| g acting on w simply pulls down the 
action from the exponent and then varies the metric leading to < T^^ >. Thus the entire 
RG evolution is governed solely by <T >. 

At fixed points of the RG flow (conformal fleld theories) the Hamiltonian vanishes, 
because the /9-functions do. However once the theory is solved and explicit expressions for 
(f)a in terms of g°'{t) and t are substituted into the Hamiltonian, the resulting function is not 
analytic at flxed points. Derivatives higher than the flrst may be singular, as the example 
of A</3^ treated in the next section shows. This is to be expected since the Hamiltonain is 
deflned in terms of the free energy density which is non-analytic at critical points. 

The analysis so far has assumed that a subtraction procedure is chosen so that the 
/9-functions only depend on k implicitly through g"'{K,) and have no explicit k, depen- 
dence. Sometimes, however, it may be convenient to use a subtraction procedure which 
results in /3-functions which have an explicit k dependence, f3°'{g,t). This can be incor- 
porated into the present framework by considering t to be like an extra coupling and 
extending the n-dimensional manifold 7W to a n + 1-dimensional manifold M.e with t 
as the extra co-ordinate. The momentum conjugate to t is (minus) the Hamiltonian, 
(/)t = dtw = —H{g,(j),t). The 2n-dimensional phase space T*{M) is now extended to a 
2n + 2-dimensional phase space T*{M.e)- By deflnition one has = 1. This is clearly 
analogous to the situation in classical mechanics where phase space is extended to include 
the energy and time as extra co-ordinates. The Hamiltonian on the extended phase space 
is 

nE{gA.t)=(3H9.t)<Pa + <Pu (61) 

and a new evolution parameter r is introduced which is ultimately identifled with t when 
the theory is solved. When the theory is solved one has He = which is just the 
Hamilton- Jacobi equation for H{g,(f),t) = /3"</>a, 

dw / dw 



^+H{,,J^,t)=0. (62) 



* The usual Legendre transform involves only one operator j{x)ip{x), composite opera- 
tors being obtained by multiple functional differentiation at the same point. But this can 
be extended to include sources for composite operators [5] . 
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On the extended phase space Hamilton's equations are supplemented by 



d(f)t 
dr 



(63) 




In other words, with t = t, ^ = ^ = dtl3°'4>a and the t dependent Hamiltonian H{g, 0, t) 

on T*(A4) is not a RG invariant when such subtraction procedures are used. Apart from 
this difference the analysis is the same as before and the evolution can be described on 
T*[M.) with explicitly t dependent /^-functions. 

Finally one might ask, what is the special ingredient of renormalisation group flow 
which allows it to be written in Hamiltonian form? After all, one cannot expect any vector 
flow to be expressible as a Hamiltonain flow. The crucial ingredient is the fact that ^ = 0. 
Thus the renormalisation group is, as the name implies, a symmetry. This leads to the 
existence of the one-form G, which reflects this symmetry, as it is constant along the RG 
flow, = since — 0. In particular, as stated at the beginning of this section, 

^ and one cannot, in general, get Hamiltonian flow if A is ignored. Note that 

restricting to surfaces of constant 9 = 0t=const gives — H, as expected for a 
Hamiltonian vector field. 
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§4 An Example - Massive \(p^ 



The construction of the previous section will now be applied to the example of mas- 
sive Xip'^. We shall use dimensional rcgularisation and minimal subraction in D = 4 — e 
dimensional Euclidean space. We start with the Lagrangian (10), with a coupling for the 
identity operator included and set to unity, 

1 1 

^0 = ^di^(pod''(po+jo(po + -mlifil + 167r^^(^^ + AqI (64) 

(the factor 167r^ is inserted to tidy up some later formulae - it is Aq which is really the 
expansion parameter in (^^ theory). It will be assumed that ttt-q > and that radiative 
corrections do not change this. A cubic term has been omitted because it can always be 
eliminated by shifting the field (fQ^(f + const, (independently of the scaling that sets 
ko = 1). 

The scalar field cpo has canonical dimension 1 — f and the couplings jo, niQ, Aq and Aq 
have canonical (mass) dimensions 3 — |, 2, e and 4 — e respectively. For simplicitly we shall 
restrict ourselves to the consideration of theories symmetric in (po —>■ —ipo and set jo = 0. 

The rcgularisation of the composite operators follows the analysis of [3] , but note that 
here Aq has the opposite sign to that reference. On dimensional grounds one expects 

Ao = K^-^(A + m^F(A,e)), (65) 

for some function F(A, e), analytic in the renormalised coupling A when e 7^ (as stated 
earlier all the renormalised couplings are rescaled to be dimensionless by multiplying by 
canonical powers of k - thus = K~^rh^ where is the usual mass with canonical 
dimension two). 

In dimensional rcgularisation the renormalised couplings are given in terms of the bare 
couplings by 

= K~'^ Z2(X, e)mQ X = K,~^z\(X,e)Xo 
A = K-4+^Ao - K-VoziiX, e)F(A, e), (66) 

where Z2 and zx are the usual renormalisation co-efficients for the mass and (p^ coupling. 
Had the linear source been retained, wave function renormalisation would also have entered 
via j = K'^^^^^zijo- 

The /^-functions are immediate. 



^' = n^ = -e + m with ^ = z-\'^ 



pm- = «^ = (_2 + S{X))m^ with 6{X) = Z2^K^ (67) 
= = (-4 + e)A + m^C(A) with C(A) = (e - 25(A))F - 
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In particular the potential factorises in minimal subtraction, U{\inn?) = m'^C,{\). These 
/3-functions can be supplemented by wave function renormalisation which may be obtained 
from the /3- function for j, 

= = (-3 + I + 7(A)) with 7(A) = z^^n'^ (68) 

and i is set to zero after 7(A) has been extracted. Of course e can be put to zero in 
all of these /3-functions since everything is finite. As is well known, the functions can be 
developed in a power series in A 

/3(A) = SA^ + o{\^) 

7(A) = ^A2 + o(A^) 

5{X) = X-\X' + o(A^) ^^^^ 
(for the A^ term in C,{X) see [14] .) From equations (66) one has 



(70) 



/_^Az^ _ m^52^ 2 m^{F(e-23)+C} ,^4-6\ 

(-eA+/3)'^ (-eA+/3)^ (-eA+/3) ^ 

22"^^^ 2m^FK^-^ 

\ K^-^ J 

Thus the renormalised operators are 

_ _ /eXz^^K'\l67v^ ^ _ /m^()z^\;^\ ifl / m^F{c - 2()) + C,}h^-' 
V-eA + /3 ) ^r^° ~ V -eX + P \ -eA + /? 

=^2"'^'v + 2m2^«^-n ^^^^ 

= K^-^l 

(strictly speaking the first of these equations should have a term on the right hand side 
involving the equations of motion for the field [3], but since these operators will only be 
used inside expectation values, this is omitted). 

In order to evaluate expectation values of the renormalised operators one has to in- 
terpret the meaning of the operators (^g and (^g- ^ functional integral approach, these 
would be the objects appearing in the exponential. When brought down by differentiation 
they should be considered as time ordered products in expectation values, but of course 
there are divergences because they consist of products of the field at the same point. These 
divergences can be treated in dimensional regularisation in the usual way, by using Wick's 
theorem for free fields and expanding expectation values perturbatively. For example for 
the operator ^\{x) Wick's theorem gives 

fM=-J£M^^h,(^^ (72) 
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where : <^o(^) • denotes normal ordering (with respect to the renormahsed mass) and 



D[.) =< ^„(.)^„(0) >= / ^j^r^y (73) 



Performing the integral using dimensional continuation yields 



Thus the expectation value is 



'<^?>=^('^Vrf-l+J)+o(A), (76) 



2 Stt V47ry V 2 

and the higher order terms can be obtained by expanding the exponential in the functional 
integral and performing the space integrals (note that <: <^o(^) •>= ^^^^ 
higher orders because the normal ordering is only defined for free fields). Including the 
identity term (71) in the expression for $^2 one finds that a pole in F cancels the pole in 
-D(O) (as it must do) and, after setting e to zero 

4 2 

< >= ^^(-1 + In™') + o(A), (76) 



where the MS scheme has been used so as to avoid —7 + ln47r terms. 
Including the next order contribution gives (after some work) 



4 4 
Km 



'^'^^ ' i„_2\ , ^(\2\\ (77) 



0^2 =< $^2 >= + Inm^ + -{{\nrr?f - Inm^) + o(A2)} 

—< $A >= «^ to all orders. 

These can be derived from the generating functional 

w{X, m^ A, t) = ^ (-^ + Inm^ + ^(1 - Inm^)^ + o(A^)) + k^A (78) 

by 0a = The calculation of the order terms in w would require extracting the 

finite part of the three loop diagram (Q)- 
The Hamiltonian is thus 

iy(^,0) = /?«(^)0„ + /?^0A 

= /?(A)0A + (-2 + 5(A))m20^2 + {m\{X) - 4A)0a. 



20 



It is not difficult to prove, from equations (67), (69) and (77), that this is a renormahsation 
group invariant, ^ = 0, to order A. In fact one finds H = = —4w when the exphcit 
solutions for the VEV's (77) are substituted into the Hamiltonian. This is a consequence 
of the fact that all couplings have been scaled to be dimensionless so that the only explicit 
K dependence in (78) is the trivial volume factor. 

The renormahsation group evolution of the vacuum expectation values is thus given 

by 



f (80) 

where the potential is U (g) = +o(A^). Using the explicit expressions for <j)a in (77) and 
the /?- functions, (67) and (69), it is a straightforward calculation to check these equations 
explicitly to order A. 

Three comments on the analysis presented here are in order. 

(i) Using dimensionless couplings will always give ^| g = Aw Thus H = —Aw. However, 
it must be stressed that in applying Hamilton's equations g°' and 0^ in H{g, (p) must be 
considered as independent vairiahles. It is only after the theory has been solved and (paiffj t) 
determined as a function of gr" and t that H can be identified with —Aw. 

(ii) The Hamiltonian, if = — < T >, is a constant of the motion (RG invariant). This 
is consisted with the observation that conservation of momentum implies that the energy- 
momentum tensor does not get renormalised (provided the subtraction procedure is com- 
patable with translational invariance [15] , as MS is). Thus the bare energy- momentum 
operator is equal to the renormalised one and < To >=< T >. (The energy-momentum 
operator here is, of course, the "improved" operator of reference [15], obtained by coupling 
the the scalar field <fo to the curvature scalar TZ of D-dimensional space before varying 
the metric and only setting TZ — afterwards.) How then can the statement that the bare 
couplings are RG invariants be reconciled with the fact that —H =< T >=< Tq >^ 0? 
One answer is that H has canonical mass dimension D despite being a RG invariant. The 
dimensionless Hamiltonian H = k~^H vanishes at high energy. Perhaps a more rigorous 
way of stating this is to observe that, as we are using dimensionless renormalised couplings 
so we could also use dimensionless bare couplings 

Ao = K^~'^Ao, rfiQ = K^rriQ, Aq = «'^Ao, (81) 

where the quantities with the tildes are the dimensionful ones. One then obtains 

dw{go,(f)o) _ a ,0 



dt 



/5o>y+ < To >= 0. (82) 



Since Pq are now non-zero (they are simply the negative of the canonical dimensions of 
the couplings) this equation is quite consistent. For a regularisation procedure which 
does not preserve translational invariance one looses the Ward identity that protects the 
energy-momentum operator from renormahsation and it is no longer necessary to have 
< To >=< T >. 
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(iii) Setting ko = 1 and then forgetting about it may at first sight seem a httle dangerous. 
After all it is clearly related to wave-function renormalisation and should run like all the 
other couplings. The fact that the operator d^ip^d^ipQ can be ignored with impunity is 
related to the equations of motion (the Schwinger-Dyson equation). There is always one 
linear combination of the basic operators $° which does not get renormalised, namely that 
corresponding to the "equations of motion" Eq{x) = (fQ (x) g^^°^y 

There is a "Ward identity" (the Schwinger-Dyson equation) which ensures that this com- 
bination of linear operators does not get renormalised, so = Eq. For this reason one 
of the operators in the original Lagrangian is always redundant and can be ignored, and 
here it is d^ipod/^ipo that has been ignored. Strictly speaking though d^ipod^ipo is not the 
same operator as —(poU (po and they should be treated separately. A complete analysis is 
given in [3] , but these complications are omitted here in the interests of clarity. 
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§5 Symmetries 



The notion of a Poisson bracket structure for the renormahsation group evohition, as 
introduced in section three, immediately raises the question of how symmetries might be 
implemented on the phase space {g'^, (f)a)- As a first, almost trivial, example of a symmetry 
consider a A'"-component scalar field (^Q) * = 1, • • • , -/V with Lagrangian 

N N 

^o^^Yl ^M^o^'^o + 5^^^°)' + 4! 5Z K'''^hvWov'o- (84) 

i=l i=l ' ijkl 

There are in general ^A^(A^+ 1)(A^ + 2)(A'^ + 3) different couplings Aq'^^ and each of these 
could renormalise differently. The different renormalisations of the various fields would 
result in N different renormalised masses as well as the renormalised couplings X^^'^'' 
(the renormalised couplings are not all independent parameters, of course, being functions 
of only ^A^(A^ + 1)(A'" + 2)(A'" + 3) + 1 bare couplings). If however the bare theory enjoys 
global SO{N) invariance all the (/Jq couplings reduce to only two, which will be denoted by 
Ao and Aq. Furthermore, if there are no anomalies, this symmetry survives at the level of 
the renormalised couplings to give only three renormalised parameters m^, A and A'. There 
is a Ward identity which demands that all the renormalised masses mf must renormalise the 
same way, i.e. mf = z^rriQ, with the same renormahsation constant z^, for all the masses 
mf. Similarly all of the renormalised Lp^ couplings (a priori ^N{N + 1)(-^ + 2)(A^ + 3) 
in number) reduce to only two, A and A'. Thus the phase space, which is in principle 
^N{N + 1){N + 2){N + 3) + 2N dimensional, is reduced to being only six dimensional by 
the symmetry. 

A less trivial example is supplied by massless QED coupled to a massless charged 
scalar field, with Lagrangian 

Co = + ii^o^D^i^o + {D^^o)\D^^o) + (85) 
where the co- variant derivatives are defined by 

Df,ipo = {d^ + zeo^oM)V'o (86) 

Dfj,(po = {d^ + ieoAofj,)(po, (87) 

and 

i^o = {^^.Ao. - a,yloM)(5'^^o - d'^A!^). (88) 

A coupling for the identity operator is not necessary since the theory is massless. Also 
no gauge fixing term is included for the moment because a perturbative analysis will not 
be used here. In order to avoid volume divergences the theory can be formulated with 
periodic boundary conditions, i.e. on a four dimensional torus T^. 

There are three independent couplings, cq, cq and Aq, but Ward identities force cq 
and Co to renormalise in the same way so that their renormalised couplings are related to 
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the bare counterparts with the same renormahsation constant (i.e. e = zcq and e = zcq 
where z is the photon wave function renormahsation constant = zA^). Thus 

which imphes that e/?*^ = e/?'^. The /3-functions are therefore not independent. In particular 

1 



which immediately suggests a change of co-ordinates from (e, e) to r = ye^+e^ and 
i& = tan-^(e/e), with 0<^< 7r/2, so that is a RG invariant, |f = 0. 
The conjugate variables to e, e and A are 



V Jrp4 

(f>e = T7 d^x < A^D^if > + Hermitian conjugate 

V Jrpi 



(91) 



where V is the volume of the torus. The integration over space is kept explicit here so as 
to mantain gauge invariance which is easily proven by integrating by parts and using the 
equations of motion. The two VEV's still have canonical mass dimension D due to the 
volume factors outside the integrals. 

Since the theory is massless the potential U vanishes and the RG equations for the 
VEV's are 

+ (ae/3^)</.e + (ae/3^)</.e + {de^^)<^X = 



dt 

+ (ae/3^)0e + {deP')cl>e + (5e/3^)0A = (92) 



dt 



However these are more elegantly expressed in the (r, -d) variables with 

0r = cos 'd + sin d 0e 

(j)^ = —r sin (f)e + r cos 0e 

and 

/3'' = cosi?/3^+sini?/3^'= 



(93) 



cos -d (94) 
= 0. 

Note that since cos?? > a positive always gives a positive /?^, in agreement with the 
statement that only a non-abelian theory can be asymptotically free. 
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The RG evolution of the VEV's in terms of these new variables is now given through 
Hamilton's equations as 



d(f)r 

'dt 

d4>i) 
dt 



+ {d^n<Pr + {d^/3^)<Px - (95) 



^ + {dxPn^r + {dxP^)(f>X = 0. 



dt 

Since P'^ vanishes the Hamiltonian 

H{r, ^, A, <Pr<px) = P'^ir, ^, X)<Pr + (3\r, ^X)<Px (96) 

is independent of 0^ and is a constant of the motion. In analogy with classical mechanics 
§ might be called an ignorable co-ordinate, but the roles of co-ordinate and momenta in 
the RG are really reversed from those of classical mechanics. It is more correct to say 
the (f)^ is an ignorable expectation value because the Hamiltonian is still a non-trivial 
function of 'd. The situation is in fact more involved here than in classical mechanics 
since the momentum dual to an ignorable co-ordinate in classical mechanics only appears 
quadratically in the Hamiltonian, whereas the 'd dependence of H in equation (96) can be 
much more complicated. The invariant i? however still plays the same role as an invariant 
in classical mechanics - since H} = it generates a symmetry on phase space via 
the Poisson bracket operation. This can be viewed as the implementation of the Ward 
identities on phase space. 

The fact that one of the expectation values can be eliminated from the phase space 
can be understood from a physical point of view in the following manner. If the gauge 
field is rescaled by T^^Ofi then the gauge coupling completely drops out of the 

matter field terms in the Lagrangian (85) and only appears in the kinetic energy term for 
the gauge fields, -t^Fq, only the ratio ■& appears in the matter field Lagrangian and this 

is a RG invariant. Defining a new variable qo = ^ and forgetting about 'd there are now 

only two expectation values to be considered, 0g = | < > and 0a- One may expect 
(pq at large energies, even in a massless theory, since the /^-function for e is positive 
[16]. 

There are further interesting aspects of gauge theories when a gauge fixing term is 
added. If a term riQ{d^AQ)^ is introduced into the Lagrangian (85) then there is another (3- 
function, jS'^ = In minimal subtraction schemes, the /3-functions for the other couplings 
are independent of the gauge fixing parameter ?], so 9^/3'' (r, "dX) = dr]P^{r, 'dX) = 0. There- 
fore the expectation value cj)^ = y Jrj.^ d^x < (d.A)'^ > evolves under RG fiow according 
to 

^ + {d,P^)cl,, = ^ cl>r,{t) = cf>,{to)e~ O'^'''^"'. (97) 

If 0^ is zero at some value of t, then it is zero at all values and 0^ is another constant 
constant of the motion. This refiects the fact that r] plays no physical role in the theory. 
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If a renormalisation prescription other than minimal subtraction is used, however, it may 
not be the case that the other couphngs have /3-functions which are independent of the 
gauge fixing parameter. However, holding to the philosophy that a change in regularisation 
prescription is just a change in co-ordinates, there must be some quantity that is a RG 
invariant, i.e. the co-ordinate transform of the dimensional regularisation co-ordinate r] - 
it may look messy in the new co-ordinates but it must exist. 

In conclusion it would seem that, just as in classical mechanics, the Hamiltonian 
framework is a very powerful one for the discussion of symmetries, which play such a 
central role in all discussions of quantum field theory. 
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§6 Conclusions 



Before summarising the main results of this paj^er a few comments will be made 
about the global topology of the space of couplings A4, as promised in the introduction. 
Consider for definiteness massless QED with the only coupling being the electron charge e 

2 

(the identity operator can be ignored). It is really ct = ^ (or l/ct) which is the important 
parameter, and a must be positive since a negative value would mean that the theory 
would be unstable, as pointed out by Dyson [17] . Thus a = is not an analytic point 
and one cannot continue, even infinitesimally, to a < 0. Similarly 1/a = cannot be 
an analytic point. The manifold Ai, the positive real line, has two boundary points both 
of which are non-analytic points of the theory. For the higher dimensional case points of 
non-analyticity are clearly also of central importance. Such points might be isolated or 
might form sub-manifolds of M. with dimension k < n. If the non-analytic points were 
to form a sub-manifold of co-dimension one (a hypersurface k = n — 1), then this would 
act as an effective boundary dA4. In any event, it is clear that any understanding of the 
global topology of A4 will be inextricably linked with an understanding of the points of 
non-analyticity. 

In summary it has been argued that the renormalisation group evolution of couplings 
and vacuum expectation values can be described as a Hamiltonian flow on the 2n dimen- 
sional phase space T*{J\4) with the Hamiltonian given by 

H{g, <t>, t) = (i\g, t)ct^a + iU{g, t) - DK)ci>^, (98) 

which can be identified with minus the vacuum expectation value of the trace of the 
energy-momentum tensor. 

The natural variables canonically conjugate to the couplings are the expectation val- 
ues, 0o = daw{g,t), where w{g^t) is the generating functional or free energy density. For 
theories with massive couplings the cosmological constant plays a central role, since its 
/^-function f3j^ = U{g,t) — DA gives rise to a potential which acts as an efi'ective force in 
the RG evolution of the VEV's 

^ = -K^dC/(^,t), (99) 

where (p = (padg"" and dU = daUdg"". 

The RG evolution of any function A{g,(l), t) can be determined from the Hamiltonian 



dA . . dA 



(100) 



In particular the RG evolution of the Hamiltonian itself is given by 

— = — = dtp-{g,t)cf>a + dtU{g,t)cl>j,, (101) 
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and H is a RG invariant if a subtraction procedure is chosen so that the /3-functions are 
independent of k. Alternatively the cosmological constant can be omitted and the RG 
evolution of the Hamiltonian h on T*{A4) is given by 

^ = ^ = dtP'-ig, t)cPa + {DU{g, t) + dtU{g, t)) , (102) 

which can only be a RG invariant if both U = and the /9-functions have no explicit k 
dependence. 

The RG equation for the /^-point Green functions is 

i=i (55) 



= (^MdbP') + K^'d.u) ^Gi'^la^ {g, cl>, t). 



The crucial ingredient that gives rise to Hamiltonian flow is the underlying symmetry of 
the renormalisation group, reflected in the fact that ^ = 0. Ward identities give rise to 
further constants of the motion which generate symmetries on phase space via the Poisson 
bracket structure. 

The RG equation for the generating functional w{g{t),t) can be interpreted as a 
Hamilton-Jacobi equation 

+ Hig,^,t) = 0. (62) 

g (^g 



dw 
'dt 



Table 1 provides a summary of the correspondence between concepts in quantum fleld 
theory or statistical mechanics and classical mechanics 

It is a pleasure to thank William Deans, Chris Ford, Bill McGlinn, Charles Nash, 
and Lochlainn O'Raifeartaigh for many useful conversations about the RG. It was Denjoe 
O'Connor and Chris Stephens who first suggested to me that a symplectic structure should 
play an important role in quantum field theory and statistical mechanics, with the VEV's 
being conjugate momenta, and that the RG equation should be viewed as a Hamilton- 
Jacobi like equation. This work would not have been possible without their prompting. 
I also wish to thank the Alexander von Humboldt foundation whose financial support 
enabled me to carry out most of the calculations of section four. 
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A Quantum Field Theory-Classical Mechanics Dictionary- 



Quantum Field Theory or 
Statistical Mechanics 

Couplings g°' (t) 
/3-functions 

Vacuum expectation values (I'ait) 
Bare couplings {goAa) 

Generating functional w{g{t), go,t) 
(Free energy density) 

Hamiltonian 

H{g,cP,t)=(3H9,t)<Pa 
Hamilton's equations 
3<^- 9H_ ^ _ _dH_ 

Potential U{g, t) = (3^ + DA 
RG flow of VEV's ^ = -dU 

at 

RG equation for w 

^ + H[git),^,t)=0 

Massless theory U = 
No explicit k dependence in /? 
Anomalous dimensions 
RG invariant {t9, H} = 
(Ward identity) 



Classical 
Mechanics 

Co-ordinates q"" (t) 

Velocities (^"'{t) 

Momenta Pa{t) 
Initial point {Qo,Pa) 

Action S{q{t),qo,t) 
(Hamilton's principal function) 

Hamiltonian 

H{q,p, t) - ^g''\q)PaPh + U {q, t) 

Hamilton's equations 
Aa _ ■ _ dH 

^ ~ dpa ~ 

Potential U{q,i) 
Newton's 2^^ ^ ^ 

Hamilton-Jacobi equation 

f + if(«W,f,*) = o 

Free particle t/ = 
Conservative system 
Pseudo-forces (Coriolis) 
Constant of motion H} = 

(Symmetry generator) 



